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Abstract 

The exotic baryon Q^{uudds) is studied with microscopic calculations in 
a quark model by using a method of antisymmetrized molecular dynam- 
ics(AMD). We predict narrow states, J"" = l/2+(/ = 0), J"" = 3/2+(/ = 0), 
and J'^ = 3/2^(/ = 1), which nearly degenerate in a low-energy region of the 
uudds system. We discuss NK decay widths and estimate them to be F < 7 
for the J'^ = {1/2+, 3/2+}, and r < 1 MeV for the J'^ = 3/2" state. 
The evidence of an exotic baryon G"*" has recently been reported by several experimental 
groups. This discovery proved the existence of the multiquark hadron, whose minimal quark 
content is uudds as given by the decay modes. The study of pentaquarks has become a hot 
subject in hadron physics. A chiral soliton model [1] predicted a narrow 0+(J'^ = 1/2+) 
state whose parity contradicts the naive quark model expectation. Theoretical studies were 
done to describe G+ by many groups [2,3]. The spin parity of G+ is not only a open problem 
but also a key property to understand the dynamics of pentaquark systems. 

In this paper we would like to clarify the mechanism of the existence of narrow pentaquark 
states. We try to extract a simple picture for the pentaquark baryon with levels, width, spin- 
parity and structure from explicit calculation. In order to achieve this goal, we study the 
pentaquark with a flux-tube model [6,7] based on strong coupling QCD, by using a AMD 
method [4,5]. 
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In the flux-tube model, the interaction energy of quarks and anti-quarks is given by the 
energy of the string-hke color-electric flux, which is proportional to the minimal length of the 
flux-tube connecting quarks and anti-quarks at long distances supplemented by perturbative 
one-gluon-exchange (OGE) interaction at short distances. For the q'^q system the flux-tube 
conflguration has an exotic topology. Fig. 1(c), in addition to an ordinary meson-baryon 
topology. Fig. 1(d). An important point is that the transition between the different flux- 
tube topologies (c) and (d) is strongly suppressed because it takes place only in higher order. 
(In 1991, Carlson and Pandharipande studied exotic hadrons in the flux-tube model [8] and 
calculated a few q'^q states with very limited quantum numbers.) 

We apply the AMD method to the flux-tube model and calculate the uudds system. The 
AMD is a variational method to solve a flnite many-fermion system. One of the advantages of 
this method is that the spatial and spin degrees of freedom for all particles are independently 
treated. This method can successfully describe various types of structure such as shell-model- 
hke structure and clustering (correlated nucleons) in nuclear physics [4,5]. With the AMD 
method we calculate all the possible spin parity states of uudds system, and analyze the 
wave function to estimate the decay widths of the obtained states with a method of reduced 
width amplitudes. 

In the present calculation, the quarks are treated as non-relativistic spin-| Fermions. We 
use a Hamiltonian as H — Ho + Hj + Hf, where Hq is the kinetic energy of the quarks, Hj 
represents the short-range OGE interaction between the quarks and Hf is the energy of the 
flux tubes. Hq and Hi are represented as follows; 



where mj (the i-th quark mass) is rUq for auor d quark and for a s quark, and Tq denotes 
the kinetic energy of the center-of-mass motion. Here, we do not take into account the mass 
difference between ud and s in the second term of Hq, for simplicity, etc is the quark-gluon 
coupling constant, and is the generator of color SU (3). In Hj, we take only the dominant 
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terms, Coulomb and color-magnetic terms, and omit other terms. 
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FIG. 1. Flux-tube configurations for confined states of qq (a), q^ (b), q'^q (c), and disconnected 
flux-tube of q^q (d). Figures (e) and (f) represent the flux tubes in the color configurations, 
and [un][d(i]s, respectively. The string potentials given by the flux tubes (b) and (c) are 
supported by Lattice QCD [9] . 

In the fiux-tube quark model [6], the confining potential is written as Hf = aLf — 
where a is the string tension, Lf is the minimum length of the flux tubes, and M° is the zero- 
point energy. M° depends on the topology of the flux tubes and is necessary to fit the qq, 
q^ and q'^q potential. In the present calculation, we adjust the M° to fit the absolute masses 
for each of 3g-baryon and pentaquark. For the meson and 3g-baryon systems, the fiux tube 
configurations are given as Fig. 1(a) and (b). For the pentaquark system, the different types 
of flux-tube configurations appear as shown in Fig.l(e),(f), and (d), which correspond to 
the states, |^(e)) — \[ud][ud]s) , |$(/)) — \[uu][dd]s) , and \^{d)) — respectively 
{[qq] is defined by color anti-triplet of qq). In the present calculation of energy variation, we 
neglect the transitions among |$(e)), |*^'(/)) and and solve 5q wave functions within the 

model space (e) or (f), which corresponds to the confined states. It is reasonable because the 
transitions are suppressed as mentioned before. In the practical calculations of the string 
potential the minimum length of the flux tubes Lf is approximated by a linear 

combination of two-body distances as Lf ^ 1(^12 + 1^23 + 1^31) for a 3g-baryon, and Lf ^ 
\ (ri2 + r34) + 1 (ri3 + ru + r23 + ^24) + 1 (rn + ris + rjg + ru) for or of the pentaquark 
systems. We note that the confinement is reasonably realized by the approximation for 
$(e,/) as follows. The fiux-tube configuration (e)(or (f)) consists of seven bonds and three 
junctions. In the limit that the length(i?) of any one bond becomes much larger than other 
bonds, the approximated {Hf) behaves as a hnear potential aR. It means that all the 
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quarks and anti-quarks are bounded by the linear potential with the tension a. Therefore, 
the approximation for $(g) or is a natural extension of the usual approximation for 3g- 
baryons. It is easily proved that the approximations are equivalent to pa 
where O = —^crJ2i<j F^F^'^ij — M°, within each of the flux-tube configurations. 

We solve the eigenstates of the Hamiltonian with a variational method in the AMD 
model space [4,5]. We take a base AMD wave function in a quark model as follows. 

$(Z) = (1 ± P)A [<i,z,<j>z, ■ ■ ■ <I>Z,X] , « e-ii^^'--^^^^)', (3) 

where 1 ± P is the parity projection operator, A is the anti-symmetrization operator, and 
the spatial part (pZi of the i-th single-particle wave function is written by a Gaussian with 
the center is a complex parameter). X is the spin-isospin-color function. For the 

pentaquark(MM(i(is) system, X is expressed as 

-X Smi,m2,m3,m4,m5 ^mim2m3m4m5\^1^2^3^i'^5) S 

<S> {\ududs) or \uudds)} <S) eabgeceh^ghf\abcef)c, (4) 

where \ududs) and \uudds) correspond to the configurations [wo?][xio?]s and [iiw][o?o?]s in Fig.l, 
respectively. Here, \a)c{a — 1,2,3) denotes the color function, and \m)s{m i) is the 
intrinsic-spin function. Since we are interested in the confined states, we adopt those model 
space for the color configurations {qq)3{qq)3q, but do not use the meson-baryon configura- 
tions {qqq)i{qq)i- The variational parameters are Z = {Zi, Z2, ■ • • , Z^} and Cmim2m3m4m5 
which specify the spatial and spin configurations. The energy variation for Z is performed 
by a frictional cooling method, and the coefficients C^im2m3m4m6 a-re determined by diago- 
nalization of Hamiltonian and norm matrices. After the energy variation, the intrinsic-spin 
and parity S"^ eigen wave function $(Z) for the lowest state is obtained for each S'^. 

In the numerical calculation, the linear and Coulomb potentials are approximated by 
seven-range Gaussians. We use the parameters, etc = 1.05, A = 0.13 fm, ruq — 0.313 GeV, 
a = 0.853 GeV/fm, and Arris = itIs — rrig = 0.2 GeV. The quark-gluon coupling constant 
etc is chosen so as to fit the N and A mass difference. The string tension a is adopted to 



adjust the excitation energy of A'"*(1520). The width parameter b is chosen to be 0.5 fm. 
By choosing Mq as Mjs = 972 MeV, the masses of N, iV*(1520) and A are fitted [10], and 
the masses of A, E and E*1385 are well reproduced with these parameters. 

Now, we apply the AMD method to the uudds system. For each spin parity, we calcu- 
late energies of the [wrfjfwcijs and states and adopt the lower one. In table.I, the 
calculated results are shown. We adjust the zero-point energy of the string potential Mq 
as Mjg = 2385 MeV to fit the absolute mass of the recently observed ©+. This M^^- for 
pentaquark system is chosen independently of M°3 for Sq'-baryon. If M°4^ = |-^°3 is assumed 
as Ref.[8],the calculated mass of the pentaquark is around 2.2 GeV, which is consistent with 
the result of Ref.[8]. 

The most striking point in the results is that the S'^ — 3/2~ and S'^ — 1/2+ states nearly 
degenerate with the S'^ = l/2~ states. The S'^ = 1/2+ correspond to J'^ = 1/2+ and 3/2+ 
with ,5 = 1/2, L = 1, and the S"" = 3/2" is J"" = 3/2-{S = 3/2, L = 0). The lowest state 

— 1/2~(S^ — l/2~, L — 0) exists just below the — 3/2" state, however, this state, as 
we discuss later, is expected to be much broader than other states. Other spin-parity states 
are much higher than these low-lying states. 

The L5'-partners, J'^ — 1/2+ and 3/2+ exactly degenerate in the present Hamiltonian 
where the spin-orbit and tensor terms are omitted. If we introduce the spin-orbit force into 
the Hamiltonian the L5'-splitting is small in the diquark structure because the effect of the 
spin-orbit force from the spin-zero diquarks is very weak as discussed in Ref . [12] . As shown 
later, since the present results show that the diquark structure is realized in the = 1/2+ 
and 3/2+ states, the L5'-splitting should not be large in the uudds system. 

Next, we analyze the spin structure of these states, and found that the J'^ — 
{1/2+, 3/2+}(S' = 1/2, L = 1) states consist of two spin-zero wrf-diquarks, while the 
J'^ — 3/2 ~ consists of a spin-zero iid-diquark and a spin-one lid-diquark. Since the spin- 
zero itd-diquark has the isospin 7 = and the spin-one itd-diquark has I — 1 because of 
the color asymmetry, the isospin of the J'" = 3/2~ state is / = 1, while the even-parity 
states — l/2+,3/2+ are / = 0. We consider that the J'^ — 1/2+ state corresponds to 
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Kinetic(//o) 


3.23 


3.22 


3.36 


3.19 


3.19 


Stvmg{HF) 


-0.67 


-0.66 


-0.55 


-0.64 


-0.64 


Coulomb 


-1.05 


-1.04 


-0.99 


-1.03 


-1.03 


Color mag. 


-0.01 


0.01 


-0.25 


0.04 


0.19 


^^Color mag. 


-0.06 


-0.01 


0.00 


0.02 


0.06 


E 


1.50 


1.53 


1.56 


1.56 


1.71 



TABLE I. Calculated masses(GeV) of the uudds system. The expectation values of the kinetic, 
string, Coulomb, color-magnetic terms, and that of the color-magnetic term in qq pairs are listed. 
The S"" = 3/2+ and S"" = 5/2+ states are higher than the S'^ = 5/2' state. 



the O+(1530) in the flavor 10-plet predicted by Diakonov et al. [1]. The odd-parity state, 
— 3/2~ has I — I, which means that this state is a member of the flavor 27-plet. We 
denote the .r = 1/2+ , 3/2+(/ = 0) by Q^, and the J" = 3/2-(/ = 1) by Qf. 

Although it is naively expected that unnatural spin parity states are much higher than 
the natural spin-parity l/2~ state, the results show the abnormal level structure of the 
{ududs) system, where the high spin state, J'" = 3/2^, and the unnatural parity states, 
J'^ — {l/2+,3/2+}, nearly degenerate just above the J'^ — l/2~ state. By analysing the 
details of these states, the abnormal level structure can be easily understood with a simple 
picture as follows. As shown in table. I, the J"^ = {1/2+, 3/2+}(S' = 1/2, L = 1) states 
have larger kinetic and string energies than the — 3/2~{S — 3/2, L = 0) and = 
1/2~{S — 1/2, L — 0) states, while the former states gain the color-magnetic interaction. It 
indicates that the degeneracy of the even-parity states with the odd-parity states is realized 
by the balance of the loss of kinetic and string energies and the gain of the color-magnetic 
interaction. In the J'^ — {l/2+,3/2+} and the 3/2~ states, the competition of the energy 
loss and gain can be simply understood from the point of view of the diquark structure as 
follows. As already mentioned by Jaffe and Wilczek [2], the relative motion between two 
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FIG. 2. q and q density distribution in the J'^ = 1/2+, 3/2+(S' = 1/2, L = 1) states of Theta+. 
The u-quark density (a), s density (b), and total quark-antiquark density (c) of the intrinsic state 
before parity projection are shown. The d-quark density is same as the u-quark density. The 
root-mean-square radius of q and q is 0.63 fm (the nucleon size is 0.5 fm). 

spin-zero diquarks must have the odd parity {L — 1) because of Pauh blocking between the 
two identical diquarks. In the J'^ — 3/2" state, one of the spin-zero wd-diquarks is broken 
to be a spin-one wd-diquark to avoid the Pauli blocking, then, the L = is allowed because 
diquarks are not identical. The L = is energetically favored in the kinetic and string 
terms, and the energy gain cancels the color-magnetic energy loss of a spin-one rtrf-diquark. 
Although we can not describe the J'^ = 1/2^ state by such a simple diquark picture, the 
competition of energy loss and gain in this state is similar to the J'^ — 3/2~. 

We remark that the existence of two spin-zero itd-diquarks in the J'^ = {l/2+,3/2+} 
states predicted by Jaffe and Wilczek [2] is actually confirmed in our ab initio calculations. 
We found that the component with two spin-zero lid-diquarks is 97% in the present J'" — 
{l/2+,3/2+} state. In Fig. 2, we show the quark and anti-quark density distributions in 
the J'" = {l/2+,3/2+} states. In the intrinsic state before parity projection, we found the 
spatial development of ud-uds clustering, which causes a parity-asymmetric shape (Fig. 2 
(c)). 

We estimate the i^iV-decay widths of these states by using a method of reduced width 
amphtudes [11]. The decay width F is estimated by the product F^, x Sfac, where F2,(a, Eth) 
is given by the penetrability of the barrier [10], and Sfada) is the probability of the decaying 
particle at the channel radius a. In the following discussion, we use the channel radius a = 1 
fm and the threshold energy Eth — 100 MeV. We here estimate the maximum values of the 
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widths, by taking into account only quark degrees of freedom. We omit the suppression of the 
transition between the confined state and the meson-baryon state due to the rearrangement 
of flux-tubes, which makes Sfac small in general. 

In case of even parity J'^ = 1/2+, 3/2+ states, the KN decay modes are the P-wave, which 
gives r°^^ a; 100 MeV fm~^. We calculate the overlap between the obtained pentaquark 
wave function and the K^n state, and evaluate the probability as Sfac = 0.034 fm~^. 
Roughly speaking, the main factors in this meson-baryon probability are the factor | from 
the color configuration, the factor | from the intrinsic spin part, and the other factor which 
arises from the spatial overlap. By using this value, the total width for K^n and K^p 
decays of the J'^ — 1/2+, 3/2+ states is estimated to be F < 7 MeV. For more quantitative 
discussions, it is important to treat the coupling with the KN continuum states, where one 
must take into account the suppression due to the rearrangement of flux-tube topologies. 

It is interesting that the KN decay width of the J'^ — Zjl" state is extremely small due 
to the £)-wave centrifugal barrier. In fact, F°=2 ~ 30 MeV fm~^ is much smaller than the 
P-wave case. Moreover, the J'^ = 3/2~{S'^ = 3/2~ ,L = 0) has no D-waye component, there- 
fore, no overlap with the KN{L — 2) states in the present calculation. Even if we introduce 
the spin-orbit or tensor forces, the KN probability (S'/ac) in the J'^ — 3/2" pentaquark state 
is expected to be minor. Consequently, the J'^ = 3/2~ state should be extremely narrow. 
If we assume the Sfac in the — 3/2~ to be half of that in the J'^ — l/2+,3/2+ states, 
the KN decay width is estimated to be F < 1 MeV. Contrary to the narrow features of the 
J'^ — 3/2~ state, in case of J'^ = l/2~, 5'-wave(L = 0) decay is allowed and this state should 
be much broader. 

In conclusion, we proposed a quark model in the framework of the AMD method, and ap- 
plied it to the uudds system. The level structure of the the uudds system and the properties 
of the low-lying states were studied. We predicted that the narrow J'^ — {l/2+,3/2+}(67=o) 
and — 3/2" (©7=1) states nearly degenerate. The widths of ©q ^-nd ©1" are estimated 
to be F < 7 MeV and F < 1 McV, respectively. Two spin-zero diquarks are found in the 
©0 , which confirms Jaffe-Wilczek picture. The origin of the novel level structure is the 5g 
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dynamics of the confined system bounded by the connected flux-tubes. We consider that 
the present results for the J"^ = {1/2+, 3/2"'"}(6j^o) states correspond to the experimental 
observation of 0"*", while the 0/=i is not observed yet. The existence of many narrow states, 
J"^ — 1/2+, 3/2+, and 3/2~, may give an light to further experimental observations. 

Concerning other pentaquarks, we give a comment on E{ddssu). The AMD calcula- 
tions indicate that the diquark structure disappears in the ddssu{l/2'^) due to the SU{3)- 
symmetry breaking in the color-magnetic interaction. As a result, the estimated width of 
the ddssu{l/2+) state is T f« 100 MeV, which is much broader than e+(l/2+). Also the 
3/2~ state is not so narrow because a 5'-wave decay channel S*(1530)7r is open. 

Finally, we would like to remind the readers that the absolute masses of the pentaquark 
in the present work are not predictions. We have an ambiguity of the zero-point energy of 
the string potential, which depends on the flux-tube topology in each of meson, 3g-baryon, 
pentaquark systems. We adjust that for the pentaquarks to reproduce the observed ©+ 
mass. To confirm the zero-point energy, experimental information for other pentaquark 
states are desired. 
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